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Abstract. In this paper we consider a generalized classical mechanics with fractional derivatives. The 
generalization is based on the time-clock randomization of momenta and coordinates taken from the con- 
ventional phase space. The fractional equations of motion are derived using the Hamiltonian formalism. 
The approach is illustrated with a simple-fractional oscillator in a free state and under an external force. 
Besides the behavior of the coupled fractional oscillators is analyzed. The natural extension of this approach 
to continuous systems is stated. The interpretation of the mechanics is discussed. 

PACS. 05.40.Fb Random walks and Levy flights - 45.20.Jj Lagrangian and Hamiltonian mechanics - 
63.50.+X Vibrational states in disordered systems 

1 Introduction The classical Hamiltonian (or Lagrangian) mechanics 

is formulated in terms of derivatives of integer order. This 



The fractional differential equations become very popu- 
lar for describing anomalous transport, diffusion-reaction 

— I ■ processes, superslow relaxation, etc [TJ[31inill] (and the ref- 
— I \ 

~^ • erences therein) . The interest is stimulated by the applica- 

>' . . . 

tions in various areas of physics, chemistry and engineer- 



technique suggests advanced methods for the analysis of 
conservative systems, while the physical world is rather 
nonconservative because of friction. The account of fric- 
tional forces in physical models increases the complexity 
in the mathematics needed to deal with them. 



ing fSllSirniSll^fTUlin]. Nevertheless, the derivation of such 

equations from some first physical principles is not an easy The fractional Hamiltonian (and Lagrangian) equa- 

matter. The fractional operator reflects intrinsic dissipa- tions of motion for the nonconservative systems were in- 

tive processes that are sufficiently complicated in nature, troduced into consideration by Riewe [T3]. His approach 

Their theoretical relationship with fractional calculus is is based on a simple look-out. If the frictional force is 

not yet ascertained in full. proportional to velocity, the functional form of a classical 
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the fractional derivative of one-half order. After applying the theory of probability. Also, we intend to investigate 
the variational technique, the obtained equation of mo- some concrete physical applications. Our paper is orga- 
tion for the nonconservative system contains a contribu- nized as follows. The appearance of the fractional deriva- 
tion of the frictional force. An extension of the fractional five in the equations of motion is conditioned on a pe- 
variational problem for constrained systems was offered culiar interaction of a physical system with environment, 
in [ni[T5] . However, as it was observed by Dreisigmeyer In Sect. [5] it is shown that the interaction is taken into 
and Young I16j. the conception conflicts with the princi- account through the temporal variable that represents a 
pie of causality. The point is that the variational principle sum of random intervals identically distributed. The prob- 
utilizes the integration by part. After this operation the ability distribution of such a sum asymptotically tends to 
left Riemann-Louville fractional derivative transforms into a stable distribution. Using this approach we derive the 
the right one which has a reversal arrow of time. Unfor- equation of motion in Sect. [3] Based on the derivation the 
tunately, the attempt to overcome this problem by means next section is devoted to the fractional oscillator. The 
of the treatment of the action as a Volterra series was effect of the fractional damping on resonance is consid- 
unsuccessful [17]. ered in Sect. [5] This feature is similar to an exponential 

A , . 1 r J i ■ • ii r i- 1 damping, but there are some differences. Further the dy- 

An alternative approach tor determining the iractional j- oi j 

,. r 1 J 1 J ii • Mo namics of coupled fractional oscillators is represented in 

equations oi motion has been developed recently m [18 . 

rrn 1-11 rxi r i- 1 i- i- Scct. 151 The forccd oscillatious of a multiple fractional sys- 

1 hey are derived by means ot the fractional normalization 

J... mi- J-,- 1. -J J 1 tem are analyzed in Sect. [71 The generalization of discrete 

condition. Ihis condition can be considered as a normal- ■' i— i o 

■ ,• I.,- r J- i -T, i i-- ■ c fractional systems on a continuous case is suggested in 

ization condition tor the distribution function m a trac- ■' °° 

^. , mi 1 1 i r ii r i- 1 1 Sect. [HI In conclusion we discuss an interpretation of the 

tional space. Ihe volume element of the fractional phase 

,. , 1 J ^. 1 , • 1 • T ii • fractional mechanics, 

space IS realized by fractional exterior derivatives. In this 

case the fractional system is described by a fractional 

power of coordinates and momenta. Therefore, such frac- 2 Operational time and subordination 

tional systems are essentially nonlinear. rpi ^- ■ a- i- i- ^ -t- • 

1 he spacetime is otten treated as a continuum. It is some- 

The purpose of the present paper is to provide a gen- thing like an ideal elastic medium. The conjecture is the 

eral method for describing the fractional systems from overarching principle of all physics. However, the concept 

some first principles. We will look at the problem from is a fiction. Nowhere does this feature of spacetime show 

another point of view. The main difference of our con- itself more clearly than at big bang and at collapse that 

sideration from above is a deep interrelation of the frac- it cannot be a continuum. The "elasticity" of spacetime 

tional temporal derivative with stable distributions from is only an approximation. There are very different view- 
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points to the space-time problem. One hopes that an ap- The important feature of time turns out to be saved, if 

propriate model for spacetime wiW be built in some sat- one chooses S(t) as a new time clock. The random pro- 

isfactory theory of quantum gravitation. Others opt for cess gives rise to the stochastic time arrow that is different 

some type of discrete structures abandoning a continuum, on the ordinary deterministic arrow |21] . 

In any case the spacetime is too rich a structure to be Although the random process is self-similar, it has 

pinned down by a single description. Perhaps, several over- neither stationary nor independent increments, and all its 

lapping and possible incompatible descriptions are need to moments are finite [2 2) . This process is non-Markovian, 

exhaust the complex variety of spacetime. but it is inverse to the continuous limit of a Markov ran- 

UsuaUy, in the ordinary mechanical description, char- dom process of temporal steps r(T), i. e. S{T{t)) — r. The 

acterizing the motion of a point particle, the time variable probability density of the process S(t) has the following 

is deterministic. Assume that the time variable represents Laplace image 

a sum of random temporal intervals T, being nonnegative i r 

P^(t,r) = ^ / (1) 

independent and identically distributed. Recall here the 

basic fact about the density of a positive stable random ^^^''e ^^notes the Bromwich path, and j = V^. 

variable. Following [H], the density g^{y) of such a vari- This probability density determines the probability to be 

able is defined by its characteristic function ^t the internal time r on the real time t [S] . Carrying out 



the change ut ^ u and denoting z = T/t"', the function 
p^(t,T) is written as t^°'Fa{z), where 



'bt 



/•oo 

/ e''^yg^{y)dy = exp{-|wr exp[-j (7ra/2) sign^]} . 

If the waiting times Ti belong to the strict domain of at- 
traction of an a-stable distribution (0 < a < 1), their sum ^ ~ 27rj 
n- V" (Ti -f r2 + • • • -t- T„ ) , n e N converges in distribution rj.^-^ -^^^^^^^ ^^^^-^^ -^^ ^^^^^^ ^^^^^1^ 3^^^ 
to a stable law. The choice of the index a in the range < Mainardi m originally attacked the probability 
a < 1 is caused by the support of the time steps T, on the density distribution for anomalous diffusion. To bend the 
nonnegative semiaxis. The continuous limit of the discrete gromwich path into the Hankel path, the function F^{z) 
counting process {Nt}t>o = max{n e N | ^Li T^ < t} is ^^^^ expanded as a Taylor series. Besides, it has the Fox' 
a first passage time. The process is conventionally denoted H-function representation ^ . Consequently, we have 



by S{t). For a fixed time it represents the first passage of 
the stochastic time evolution above that time level. The 



{l~a,a)\ ^ {-zf 



(0,1) ; ^^fc!r(l-a(l + fc))' 



random process S'(i) is non-decreasing and depends on the where r[x) is the usual gamma function. The function 
true time [20] . As is well known from the everyday expe- Fa{z) is entire and non- negative in z > 0. It vanishes 
rience, time is always running from the past to the future, exponentially for large positive z. Taking into account the 
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normalization relation Fa{z) dz = 1, it will readily handling problems in the theory of small oscillations and 

be seen that the fmiction p^{t,T) is really a probability in studying the dynamics of a rigid body. Hamiltonian me- 

density. chanics contains Lagrangian mechanics as a special case. 

In the theory of anomalous diffusion the random pro- In addition the Hamiltonian formalism permits ones to 

cess S{t) is used for the subordination of other random solve some problems (as the attraction by two stationary 

processes. Recall that a subordinated process Y{U{t)) is centers and the determination of geodesies on the triaxial 

obtained by randomizing the time clock of a random pro- ellipsoid) which do not yield solutions by other means. It 

cess Y{t) using a random process U{t) called the direct- is expected that this route may be followed in the analysis 

ing process. The latter process is also often referred to of fractional systems. 

as the randomized time or operational time |1S]. These Let a Hamiltonian system evolution depend on oper- 

concepts are helpful in formulating the continuous time ational time t. The corresponding equation of motion is 

random walk approach, deriving the equation of subdifFu- written as 
sion, finding its solutions and moments, obtaining a sep- 

aration ansatz and in many others. The various features '^'^ ^'^ '^1 

were exhaustively analyzed in literature (see, e. g., dSUHl Consider a dynamical system for which the momentum 

[^[^[5011511 13211331 134)'). The mathematically elegant sur- and the coordinate satisfy the relations 
vey is represented in the paper [22) that especially is rec- 
ommended for a deep insight into the interrelation be- 
tween fractional calculus and the theory of probability. 
Nevertheless, it should be noticed that the exact physical 
nature of the power-law waiting-time distributions in the 
subordination scheme is unclear. The problem requires a 

more rigorous examination elsewhere. Pa(.s) = e Pa{t) dt ^ s p{s ) , 

POO 

Jo 

3 Equations of motion 

where a is in the range of < a < 1. Since the values 

Many important problems from classical mechanics can dn/dp and dn/dq depend on operational time, we 



Pa{t)^ / p^{t,T)p{T)dT, 

Jo 

/•oo 

qc.it) = / /(t,T)g(r)dr. 



Their Laplace transform with respect to time has a simple 
algebraic form 



be solved using the Newtonian formalism. Among them assume that 
is the problem of motion in a central force field which is 
of fundamental importance in celestial mechanics. On the ^^P' 
other hand the Lagrangian formalism is more suited for 



dn 



- = P (t, V—dT, 
Jo dp 

P^{t,T)—dT. 
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In this case Eqs.® become fractional: the Hamiltonian Q is not an explicit function of time, 

d°^<la dT-La d'^Pa dHa non-intcgcr values a the system is nonconservative 

because of the fractional derivative of momentum. Then 
Here we use the so-called Caputo derivative [35l|3S[37] de- Hamiltonian equations for the fractional oscillator are 

fi'^^'i written as 

dt" U qa- -7— - -TT— - Pa , [d) 

' ^ ■ dr, n - 1 < a < n, fta _ _ oHg _ 2, 



r{n~a) {t-T)"+^-" D p^ = ^^-^^-u q^. (6) 

where = £>" xfi) means the n-derivative oi xlt), ^ r 

It follows from this that 

and 

J^yit) = ^^j\ir)it-rr^dr + - V = or D^,. + . V = 0. (7) 

is the fractional integral. In fact, the power-law waiting Each of the equations has two independent solutions. Suf- 

times imply fractional derivatives in time. The reader who ^ce it to solve one of these equations, for example, that 

desires further background information on fractional cal- determines the coordinate: 
cuius will do well to consult the excellent books [551139] . 



When a = 1, the generalized equations ([3]) transform 
into the ordinary, namely Hamiltonian equations being 



where A, B are constants, and 



oo ^ 

well known from classical mechanics. Thus this point of t7 i „\ _ ^ ,, ,, ^ n 

view embraces a wide circle of physical tasks. In next sec- 

is the two-parameter Mittag-Leffler function. Then the 

tions we will elaborate some of them. 

fractional oscillator momentum is expressed in terms of 
the fractional derivative of the coordinate 



4 Fractional oscillator 



Pa(t) = 7n]D°'qa{t) , 



One of the simplest physical models (but nontrivial) sup- 

ported by the above-mentioned method is a fractional os- ^'^^^^ t^*' generalized mass gD]. In this connection 



cillator. Its generalized Hamiltonian takes the form 



it is relevant to remark that 



Ha = ipl+LJ^ql)/2, (4) £;2a,l(-i ") - -t"£^2a,l+a(-t 

where a; is the circular frequency, and Pa the displace- 
ment and the momentum respectively. The value describes The phase portrait is a close curve only for the harmonic 
the total energy of this dynamical system [40 . Although oscillator {2a — 2), but in our case this is a spiral. The 
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total energy decreases. An example of the phase plane di- following equation 
agrani for the fractional oscillator is represented in Fig. [TJ ujq /"* ,/^a-l 

The intrinsic dissipation in the fractional oscihator is cau- 



1 f* 







{t - t')"-^ F{t') dt' , (9) 



sed by the following reason. The fractional oscillator may ^(o^) Jo 

be considered as an ensemble average of harmonic oscilla- where it should be kept 1 < a < 2, and F is the exter- 

tors gj. The oscillators differ shghtly from each other in ^al force. The dynamic response of the driven fractional 

frequency because of the subordination. Therefore, even if oscdlator was mvestigated m gO]: 

they start in phase, after a time the oscillators wiU be al- x{t) ^ [ F{t') {t - t')"''^ Ea,a{-uJait - t')") dt' . (10) 

Jo 

located uniformly up to the clock-face. Although each os- 



cillator is conservative, the system of such oscillators with 
the dynamics like a fractional oscillator demonstrates a 
dissipative process stochastic by nature. 

As is shown in |37lH2] . the fractional oscillations ex- 
hibit a finite number of damped oscillations with an al- 
gebraic decay (Fig. [2]) . The fact rests on that such oscil- 
lations may be decomposed into two parts. One of them 
demonstrates asymptotically an algebraic (monotonic) de- 
cay, and another term represents an exponentially damped 
harmonic oscillation. Owing to the second term decreas- 
ing faster than this happens for the term with an algebraic 

,(i) ^ J_ / e- %/"^^+/"^'^\, . (11) 
decay, the fractional oscillations possess a finite number of Sttj' (s + uj j (s" -I- Wg ) 

zeros. Our notation named as a circular frequency is The Bromwich integral ([TT]) can be evaluated in terms of 

enough acceptable, since the value characterizes the num- the theory of complex variables. Some particular exam- 

ber of oscillations on a given temporal interval. ples of driving functions in this context were considered 

in |40| . However, the case under interest was studied in 



This allows us to define the response for any desired forc- 
ing function F{t). The "free" and "forced" oscillations of 
such a fractional oscillator depend on the index a. How- 
ever, in the first case the damping is characterized only 
by the "natural frequency" ujq, whereas the damping in 
the case of "forced" oscillations depends also on the driv- 
ing frequency w. Each of these cases has a characteristic 
algebraic tail associated with damping |43) . 

If F{t) is periodic, namely F{t) = Foe-'"*, then the 
solution of Eq. Q is determined by taking the inverse 
Laplace transform 



5 Effect of fractional damping on resonance 



Now we consider a behavior of the fractional oscillator 
under an external force. Following the initial conditions 

x(0) = and i(0) — 0, this model is described by the ^' ~ ~J\a) Jg 



PT| . Consider only the harmonically forced oscillation. If 
one waits for a long enough time, the normal mode in this 
system is damped. After the substitution xqb^'^* for x{t) 
in Eq.® we obtain 



e^^' = f\t - t'r~' xo e^""*' dt' + 
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Fig. 1. Phase plane diagram of the fractional oscillator with the initial conditions qa(0) — 1, D°'qa{0) — 0, where a = 1. 



1 



{t^t')"-^ Foe^'^'' dt' . (12) 



^(«) Jo 

It is convenient to change the variable u)[t — t') = Q in 
the integrand. Next we can divide out exp(jwt) from each 
side of (|12p and direct t to infinity. The procedure permits 
one to extract the contribution of steady-state oscillations. 
Then Eq. ^ gives 



^" [u^ + exp(j^a/2)] " 
The forced solution pexp{jujt + 9) is written as 

^0 

[wg" + a;2" + 2w^a;" cos(7ra/2)] 



(13) 



P 



tan 9 = — 



w"sin(7ra/2) 
U!^ -fw" cos(7ra/2) 



Denote uj/ujq — z. Then we have 



„2, ,2a 
P ^0 

Fo 

tan 9 



1 



[z2" + 2z" cos(7rQ;/2) + 1] ' 
z" sin(7rQ;/2) 



1 + cos(7rQ;/2) 



(14) 



It should be noticed that the maximum of p^uj"^" /Fq is 
not attained for lo — ujq. It is shifted to the origin of 
coordinates on Fig. [3] with decaying a. To differentiate 
P^ujq'^ /Fq with respect to z, we get the maximum Zmax = 



\/cos[7r(2 — a)/2\. For a — 2 the damping in this oscillator 
vanishes. If then both frequencies coincide u = u^, the 
amplitude of the oscillator tends to infinity. 



8 
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Fig. 2. Fractional oscillations are characterized by a finite number of zeros {a = 1.65). 



6 Coupled fractional oscillators 



From the theory of vibrations it is well known, much of 
physical systems permits a description in the form of free 
harmonic oscillators. However, such a representation is too 
idealized, as in many practical cases these systems are not 
usually isolated. Instead, they interact with environment, 
or with other oscillators. Therefore, the study of the dy- 
namics of driven and coupled oscillators is of major im- 
portance. 

Here we intend to provide a similar analysis for frac- 
tional systems. Consider two identical fractional oscilla- 
tors mutually coupled. For 1 < a < 2 the dynamics of 



this system is given by 

D"xi+Ljlxi+K^{xi-X2) =0, (15) 
D°'X2+Ujf,X2+K^{x2~Xi) =0, (16) 

where the two oscillators are labelled by 1 and 2, respec- 
tively, and is the measure of the coupling, luq the circu- 
lar frequency. Let us introduce new variables ui — xi+X2 
and U2 ~ xi ~ X2. In such coordinates the system of equa- 
tions ((T5|) , (fT6|) transforms into the equations of two inde- 
pendent fractional oscillators with the frequencies ujo and 



^/u!q + 2k'^. If U2 — 0, then xi — X2 so that both oscilla- 
tors move in phase with the frequency ujq . In this case the 
coupling between the oscillators has no influence on their 
motion. If lii — or xi ~ ~X2 that is the same, the os- 
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Fig. 3. Resonance curves and phase shift of the fractional oscillator forced by an harmonic oscillation with the frequency 
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dilators evolve in antiphase with the increased frequency which are simple fractional oscillators, is a major result in 

a/wq + in force of the measure of the coupling. the physics of coupled fractional oscillators. This decom- 

It is often found that the physical systems of interest position is not restricted to just two coupled oscillators, 

consist of two and more oscillators which interact weekly An analogous decomposition holds for an arbitrary large 

among them. When the coupling is weak (k ^ 1), the frac- number of coupled fractional oscillators, 
tional oscillators transfer their energy from ones to others 

and vice versa. The effect depends on the magnitude of the 7 Forced OSCillations of a multiple fractional 

index a characterizing a strength of dissipation as well as 

system 

on the magnitude of the parameter k. Nevertheless, due 

^ ^, ,. . ^. „ „ ,1 r i- 1 -11 i- -11 Let a fractional oscillator have two degree of freedom. The 
to the dissipation, hnally the fractional oscillations will " 

, . i-i 1 harmonic force with the frequency p acts on the system 

decrease m amplitude. 1 j j 

, .,, , , .... 11 r,^, n,i . coordinates x and v. Then the equations of motion are 

Let two oscillators rest initially. 1 hen one of them gets 

i?"/^xi(0) — Bq. Following Sect. HI the time evolution of 

{AD°'x + HD^y + ax + hy = X cos pi , 
(19) 
HD"x + BD°'y + hx + by = Y cos pt 

B r 

xi{t) = — woi"^^£^Q.i+a/2(-'^o^") + with 1 < a < 2. We shah seek solutions to Eqs. (HH) in 



2 

+ t"'^E^^ i+a/2(-(^o + 2«')i")l , (17) 



the form 



B r 

X2{t) ^ -Z- t^Qt°'/'^Ea,l+a/2{~t^lt°') - X = (3cOSpt, y = J COS pt . 



2 

- i"/'i?a, l+a/2(-(c^0 + 2^)^)] . (18) 

Fig. m shows a superposition of two fractional oscillations. 
The oscillations observable are complicated, but they can 
be likened to damping beats. In the force of this damping 
the minimum of xi does not coincide with the maximum 
of X2 ■ The phase shift between them is determined by the 



The treatment of dynamical equations can be greatly sim- 
plified by using the representation in complex numbers. 
With this object in view, replace cos pt by e-"'* in (fT9|) . 
Next we substitute x = /3e^^* and y = je^^* in such a 
system of equations. Pick out e-'^* in each term of the 
equations. The fractional derivative of e-"'* gives 



value a, i. e. by the level of dissipation in this system. D"^e'^* = J^""_D^e-"'* — —p^J^ "gip* 



Moreover, the algebraic relaxation term of the fractional _ P ^ / _m-i -jr j 

oscillations deforms the position of the maxima of xi and 



r{m) Jo 

where m = 2 — a is constant. Making t — >■ cx) in the above 
integral, we apply the table integral |44) : 



X2- Distinctly this can be seen in Fig. I^b . 

The decomposition of the dynamics of coupled frac- 
tional oscillators in a superposition of normal modes. 



12 
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(20) 



for < Re^ < 1. This technique transforms Eqs. (TT^ into 

the algebraic expressions 

(a + Ap^eJ''^"/^)^ + {h + iJp"e^''^"/2)7 = X , 

{h + Hp°'e^''"/'^)I3 + {b + Bp"e^''°'/^)j = Y . 
The steady-state solutions of Eqs. ([T5| can be derived from 
((20|) by taking the "real" part of /3 and 7, i. e. by projecting 
the motion onto the real axis of x and y. From Eas. ipU)) 
at once it follows that 

X h + Hp"e^''°'/'^ 
^ h + Hp^e^""^^^ Y 

where 



^4 



1 



7 = 



(21) 



A = 



a + Ap^eJ'^"/^ h + Hp°'e^-""^^ 
h + Hp"e^''"^^ b + Bp"e^''"^^ 
is the determinant. Thus the general conclusion is the fol- 
lowing. When a periodic force of simple-harmonic type 
acts on any part of the system, every part executes a 
simple-harmonic vibration of the same period, but the 
amplitude will be different in different parts. When the 
period of the forced vibration nearly coincides with one 
of the free modes, an abnormal amplitude of forced vibra- 
tion will be in general result, owing to the smallness of the 
denominator on the formulae (j2ip . 
If X 7^ 0, r = 0, then 

X{h + Hp°'e^-""^^) 



1 



(22) 



Imagine now a second case of forced vibration in which 
X = 0, r 7^ 0. This yields 

^_ r(/i + iJp"e^''^"/2) 



Comparing with (0^ . we see that 

7:X = /3:y. (23) 

The result concerns a remarkable theorem of reciprocity, 
first proved for the theory of aerials by Helmholtz, and 
afterwards greatly extended by Rayleigh. Their interpre- 
tation is most easily expressed when the "forces" X, Y 
are of the same character in such a way that we may put 
X = Y and obtain /? = 7. The coupled fractional os- 
cillators forced by harmonic oscillations support well the 
theorem. Indeed, it is clear in view of their equations of 
motion being linear. 

The resonance curves are represented in Fig. [5] Here 
we consider the case, when X ^Q,Y = {]. The curves de- 
niostrate the following interesting effects: 1) if the exterior 
force frequency Q coincides with one of normal frequencies 
of this system, the amplitudes of both oscillators increase; 
2) if the frequency of the exterior force, acting on the first 



oscillator, coincides with cos(7r(2 — a)/2)} > 0, 

then the first oscillator is rest (/3 = 0). The latter phe- 
nomenon is called the dynamic damping often used for 
the breaking of unfavourable oscillations . For the ex- 



terior force frequency ^h/{H cos{tt{2 — q;)/2)} > the 
second oscillator will be rest. 

The extension of the method to the general case of lin- 
ear fractional oscillations of a multiple system is obvious 
so that the result may be stated formally. In such a system 
with k degrees of freedom there are in general k distinct 
"normal modes" of free fractional vibrations about a con- 
figuration of stable equilibrium, the frequencies of which 
are given by a symmetrical determinantal equation of the 
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Fig. 5. Resonance curves of the fractional oscillator forced with two degree of freedom: X — 3, Y — 0, a = 1.99, H = h = 
1,A = B — 2,a = b = 3. The dynamical damping is detectable, when P — and 7 7^ 0. See text for more details. 

fc-th order in p°'e^'""^'^, analogous to ^U^. In each of these ing chain as a limiting form of the problem where a large 
modes the system oscillates exactly as if it had only one number of equal and equidistant particles are attached to 
degree of freedom. a tense string whose own mass is neglected. The general 

principles may be applied to fractional systems too. In this 
section we shall be concerned with fractional systems for 
which the number of degrees of freedom becomes infinite. 



8 Transition to continuous systems 



At least mathematically, it is sometimes possible to pass The lattice represents the most illustrative example 

from the study of oscillatory systems of finite freedom which is naturally called by the ordered structure of oscil- 

to continuous systems by a sort of limiting process. So lators. To start with an one-dimensional chain of identical 

D. Bernoulli (1732) considered the vibrations of a hang- fractional oscillators, the characteristic spatial period of 
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wave motion in this chain is assumed to be too more than circadian rhythms and beats of the heart, electromagnetic 
the mesh dimension. If the oscillators interact with the waves and modes of oscillation of the atom nucleus are 
nearest neighbours, they are described by the equations examples of this kind of phenomena among many oth- 
of type ers. Electrical and mechanical oscillators are everyday con- 

J, stituents in the world of engineering. Usually an oscilla- 

+ = M{(j)n^i + (t>n+l - 20„) , (24) 

tor is something that behaves cyclically, changing in some 
where M is constant, and 1 < a < 2. Proceeding from ^ay, but eventually getting back to where it started again, 
the discrete system to a continuous one, we arrive at the The fractional calculus extends our representation about 
equation in partial derivatives oscillatory phenomena. The temporal evolution of frac- 

-2 9^4' 9 tional oscillator models occupies an intermediate place 

a^-- a^ + -o0 = o, (25) 

between exponential relaxation and pure harmonic oscilla- 
where = Ma\ and a is the mesh size. If = 0, the r^^^ fractional osciUator clearly demonstrates a pure 

equation m is simplified so that its solution is algebraic decay together with an exponentially damped 



Cl>{x,t)=l^ pS{t,T){MvT-x) + f2{vT-x)}dT, (26) 

where the functions /i, /2 are arbitrary. The two terms 
in (1261) admit of a simple interpretation. They represent 
damped-in-time waveforms traveling in the direction of 
x-positive and x-negative. 

In this connection it should be necessarily mentioned 
that for < a < 1 Eq. ([25]) will describe an anomalous 
diffusion with a potential. Then the function (j){x, t) be- 
comes a probability density. Thus the value has quite an 
other interpretation [46^ in comparison with the above- 
mentioned case. 

9 Interpretation and discussion 



harmonic motion. 

The modification of the conventional representation of 
the Hamilton equations is conditioned on a random in- 
teraction of the subsystems with environment. Each sub- 
system is governed by its own internal clock. Although 
its dynamics is described by the ordinary Hamiltonian 
equations, their coordinates and momenta depend on the 
operational time. The passage from the operational time 
to the physical time through the averaging procedure ac- 
counts for the interaction of the subsystem with environ- 
ment. Consequently, the whole of subsystems behaves as 
a fractional system. 



Nature is permeated with oscillatory phenomena. Pulsat- 
ing stars and earthquakes, oscillating chemical reactions The author thanks D. Dreisigmeyer for fruitful discus- 
and long term variations of the Earth's magnetic fields, sions. 
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